Strong-coupling behaviour in discrete 
Kardar-Parisi-Zhang equations 



T. J. Newman and A. J. Bray 
Department of Physics, 
University of Manchester, 
Manchester, M13 9PL, UK 

Abstract 

We present a systematic discretization scheme for the Kardar-Parisi-Zhang 
(KPZ) equation, which correctly captures the strong-coupling properties of 
the continuum model. In particular we show that the scheme contains no 
finite-time singularities in contrast to conventional schemes. The implications 
of these results to i) previous numerical integration of the KPZ equation, and 
ii) the non-trivial diversity of universality classes for discrete models of 'KPZ- 
type' are examined. The new scheme makes the strong-coupling physics of 
the KPZ equation more transparent than the original continuum version and 
allows the possibility of building new continuum models which may be easier 
to analyse in the strong-coupling regime. 
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I. INTRODUCTION 



Within the field of non-equilibrium interface growth, much attention has been given to 
a continuum Langevin description first suggested by Kardar, Parisi and Zhang (KPZ) [|T],Q . 
The KPZ equation for the interface profile h(x, t) takes the form 

d t h = vV 2 h + \{Whf + T], (1) 

where T](pc,t) is a gaussian white noise with zero mean. One is interested in asymptotic 
scaling behaviour of this model. In particular the steady-state two-point correlation function 
C(|r — r'\,t — t') = ((h(r,t) — h(r',t')) 2 ) is expected to assume a scaling form 

C(R,r)~I?Xf(R*/T), (2) 

where % an d z are the roughness and dynamic exponents respectively. There exists a scaling 
relation between these exponents of the form x+<2 = 2. The determination of these exponents 
along with the scaling function f(y) are the prime goals in this field. For substrate dimension 
d < 2 the exponents take on non-trivial values for all A 7^ 0. For d > 2 these exists a phase 
transition between smooth and rough phases. For small A the interface is smooth (x = 0) 
- this behaviour is described by the so-called weak-coupling fixed-point. For large A the 
interface is asymptotically rough and x takes on a non-trivial ^-dependent value - this is 
the physics of the strong-coupling (SC) fixed point. 

It is convenient for us to mention at this point the interesting connection between the 
KPZ equation and the physics of directed polymers in a random medium. This may be 
realised by making the Hopf-Cole transformation /i(x, t) = {v/X) lnw(x, t) which yields the 
equation 

d t w = uV 2 w + (A/V) w i]. (3) 

The field w is interpreted as the restricted partition function of the directed polymer (which 
is more clear if one rewrites the above equation as a Feynman path integral.) Although the 
equation is linear, the noise is now multiplicative which forbids any simple analytic progress. 
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In fact there are no systematic methods available to study the SC behaviour of the KPZ 
equation analytically. The most useful technique to date is a mode-coupling scheme which 
certainly gives non-trivial behaviour in higher dimensions, but suffers from the draw-back 
that it is very difficult to correct systematically. (In particular the mode-coupling scheme 
indicates d = 4 as an upper critical dimension above which the exponents take on mean-field 
values; i.e. x — 0.) 

In such a situation, numerical methods take on a supreme significance. In attempting 
to understand the KPZ model, two main numerical methods have been employed. The 
first is that of simulating microscopically motivated models which are believed on physical 
grounds to lie within the KPZ universality class. Such models include Eden growth [|J, 
polynuclear growth models 0, the hypercube stacking model || and the Kim-Kosterlitz 
model J7j. The dynamical exponents measured for these different models all lie within 
some small range, although numerical precision is not good enough (especially in higher 
dimensions) to give a definitive answer as to whether these models share the same universal 
features. The more subtle question of whether their universal features are actually those 
of the KPZ equation, relies upon measuring the exponents (and in principle the scaling 
functions) for the KPZ equation itself. The lack of analytic techniques means that such 
measurements may only be made by direct numerical integration of a discretized form of the 
KPZ equation. (Alternatively one may perform numerical studies of the directed polymer 
analogue via transfer matrix methods.) Such a numerical procedure has been carried out 
by several groups over the years pHlOfl with varying degrees of success. A common feature 



of these attempts is that numerical integration becomes increasingly difficult as the value 
of the coupling constant is raised. (The instability for large coupling was also noted in ref. 



IIJ.) This is unfortunate as ideally one would like to be deep inside the SC region, which 
for higher dimensions implies that the dimensionless coupling constant is very large. In 
this paper we shall show that the choice of spatial discretization of the KPZ equation is 
crucial in order to ensure stability under integration, and more importantly to ensure that 
the discrete model is still within the universality class of the original continuum model. The 



discretization schemes used in previous studies unfortunately fail on both these counts and 
the results so obtained must be viewed with scepticism. 

The remainder of this paper takes the following form. In the next section we study 
in some detail the conventional discretization schemes used to mimic the continuum KPZ 
equation. We shall show that they contain defects such as instabilities and ghost fixed points 
which are absent in the continuum. In section 3 we motivate a new discretization scheme 
which is guaranteed to be stable and which also contains no unphysical fixed points. This 
model is the ideal basis for a fresh numerical investigation of the KPZ equation. In section 4 
we look at the SC properties of this discrete model and find that it is physically more trans- 
parent than its continuum counterpart. Using this new insight we make some speculative 
propositions concerning new continuum models which are intrinsically SC versions of KPZ. 
The hope is that these models are more analytically susceptible than the KPZ equation. We 
end with section 5, in which we discuss the nature of universality in the space of discrete 
'KPZ-type' models, along with giving our conclusions. 

II. CONVENTIONAL DISCRETIZATION SCHEMES 

In order to perform numerical integration of the KPZ equation it is necessary to discretize 
both space and time. In this paper we shall be exclusively concerned with the delicate nature 
of the spatial discretization. Temporal discretization is a less challenging problem, and 
there exist many techniques (such as predictor-corrector methods) which help to stabilize 
the integration. So henceforth we will keep time continuous. Now the standard spatial 
discretization for the KPZ equation is performed as follows. One replaces the function 
h(x,t) by a set of fields hi(t) where the subscript is a lattice index. The lattice is taken to 
be hypercubic with spacing a. One then rewrites Eq.(|T]) as 

dhi(t)/dt = (v/a 2 ) - hi(t)) + (X/2a 2 ) Y^hjit) - Ht)Y + Vi(t), (4) 

3 3 

where the sums are over nearest neighbour sites. 
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There is a minor subtlety here which we shall address immediately before continuing. In 
writing the discrete form of (Vh) 2 there are two obvious choices. The one favoured in past 



At this level of the discussion we state that the second form is preferable, as the first form 
has the undesirable lattice feature that it vanishes when the height at site % is a symmetric 
peak. Studying surface morphologies generated by the first form shows that such sites get 
'left behind' by their neighbours, as this site receives no driving force from the nonlinearity. 
This point is academic in the light of the remaining discussion in this section. 

We shall take two approaches in order to illustrate the problems with the above dis- 
cretization scheme. Both approaches are concerned with the deterministic version of the 
problem, and both yield the result that even in the absence of noise the above discretization 
scheme is unstable, which is entirely unphysical as the deterministic continuum model is 
completely stable and has a single large time solution of h = const, (for bounded intial 
data.) In the first approach we compare the numerically integrated solution of Eq. @ with 
the exact solution of the deterministic continuum equation, for the simple initial condition 
of a top hat. In the second approach we exactly solve the discrete equation for some simple 
cases and show explicitly that the solutions are unstable. The inherent instability becomes 
progressively worse as the coupling A is increased in value. 

So, beginning with the first approach we consider the deterministic version of the KPZ 
equation. An exact solution of this equation may be obtained via the Hopf-Cole transfor- 
mation and takes the form 





(5) 



whilst the alternative is 




(6) 





5 



where g(x, t) is the heat kernel. 

We consider d — 1 for simplicity. Consider the simple initial condition of a top hat of 
height H and width 2b centred at the origin. The above solution then simplifies to 



M*,t) = ^ln{l 1 ^ XH ^ 



A 2 



. / b — x \ ( b + x 

erf 7- , N1/n + erf 



(Art) 1 / 2 } W4i4)V 



2 



(8) 



where erf (z) is the error function [I2|. It is easy to convince oneself that h(x,t) is a mono- 



tonically decreasing function of time for \x\ < b. In other words, the block decays away to 
zero, for any size of the coupling constant A. If one now uses Eq.(|]) to evolve the top hat 
function one finds that this smooth decay only occurs when A is small. In fact the value 
of the coupling where unphysical behaviour sets in is of order A c = v/H. In fig.l we plot 
the evolution of the height at the centre of the top hat and at one edge as a function of 
time, for A = 10A C - curves corresponding to the exact solution given in Eq.(|8]), and direct 
numerical integration of the discrete scheme of Eq.(f|) are shown. One sees that the discrete 
equation fails completely to capture the correct evolution. This effect cannot be corrected 
by reducing the time step. It is a consequence of the spatial discretization and would occur 
even for continuous time. 

This result is very disappointing as A c demarcates the transition between weak and 
strong coupling behaviour for this simple problem. From the exact solution one sees that 
interesting and revealing KPZ physics is present. For A A c the nonlinearity is essentially 
unimportant and the block diffuses away. More precisely, the centre of the block decays as 

This power-law decay also holds for A ^> A c but only for extremely large times. When 
the coupling exceeds its critical value, a new time scale t* emerges, which takes the form 
t* = jixv) exp(2A/A c ). The block diffuses away only for t > t*. For t < t* the centre of 
the block is essentially frozen, due to the nonlinearity pushing out the sides of the block. In 
fact the decay of the centre in this time regime is of the form 

h(0,t) ~ H- (v/2\)\n(irvt/b 2 ). (10) 



In Appendix A we discuss in more detail the relevance of this result to the SC physics of 
the noisy KPZ equation. 

The discrete equation fails for the top hat intial condition due to the large height devi- 
ation which occurs over a few lattice sites. We have tested the discrete equation on smooth 
initial conditions such as an inverted parabola and a gaussian, and it reproduces the correct 
evolution. In the stochastic version of the model, the noise is constantly producing discon- 
tinuities into the interface, so any discrete scheme must be stable in the presence of large, 
local height deviations. 

So to summarize the result of the first approach, we see that even for a simple scenario - 
that of the deterministic evolution of a block profile - the discrete equation fails to describe 
the strong-coupling physics. 

We now turn to the second approach - that of solving the discrete equation (|]) explicitly 
(again in the deterministic version) for some simple situations. As a warm-up we consider 
the ostensibly artificial problem of just three sites on a ring. We are led to consider the 
coupled equations 

dh/dt = h 3 + h 2 ~ 2hi + (X/2)[(h 1 - h 3 ) 2 + (h 2 - h) 2 ] 

dh 2 /dt = h 1 + h 3 - 2h 2 + (A/2) [{h 2 -h 1 ) 2 + {h 3 - h 2 ) 2 } (11) 

dh 3 /dt = h 2 + h 1 - 2h 3 + (\/2)[(h 3 - h 2 ) 2 + (h x - h 3 ) 2 } 

where we have set v — 1 = a. Eliminating the soft mode (the mean height) we have two 
independent modes which we combine as U\ = (A/6)(/ii — h 2 ) and u 2 = (X/6)(h 2 — h 3 ). The 
above equations now read (after rescaling t — > t/3) 

dui/dt = u\( — 1 + u\ + 2u 2 ) 

du 2 j ' dt = u 2 ( — l — u 2 — 2u\) (12) 

The alternative discretization of (Vh) 2 yields the same equations as above, but with A — > 
— 2A. Given this is a deterministic system, we expect on the basis of the deterministic KPZ 
equation for the final state to be that of all three heights equal, independent of the initial 

7 



condition, i.e. Ui(t) — > 0. That this is not the case for this discrete model may be seen 
by finding the fixed points of Eqs. (|T2|), along with their associated stability. There are 
actually four fixed points, only one of which is the physically reasonable one. The remaining 
three fixed points correspond to one particular configuration, allowing for the threefold 
translational degeneracy. In fig. 2 we show the dynamic flows of the above equations and 
sketch the approximate boundaries of stability. The insets shows the configurations relating 
to the various fixed points. One sees that the model is only driven to the correct fixed point 
for initial conditions which lie within the boundaries of stability. The linear dimensions of 
the region (in terms of the original hi variables) are of order 1/A. This indicates once again 
that as A is increased the stability of the discrete equations decreases. 

Although this problem seems far from any interesting situation, this is not the case. Since 
we have enforced periodic boundary conditions, the above result shows that any periodic 
chain whose length is a multiple of three will be unstable since the chain only has to partition 
itself into groups of three for the unstable fixed points found above to become active. Of 
course, as the number of particles increases, there will be an increasingly large number of 
fixed points corresponding to different modes of the system. That there exists at least one 
such unstable mode (as we have demonstrated above) is sufficient for our purposes. This 
result states that the discrete scheme given in Eq.(f§) is intrinsically unstable due to the 
existence of ghost (i.e. unphysical) fixed points which separate the evolution of the system 
from its true asymptotic state. 

To illustrate this point further, we briefly consider one more model solution of the discrete 
equations, this time for an arbitrarily large number of sites. Consider Eq.([|) for N height 
variables {hi(t)} coupled together via the scheme given in Eq.fl5|). (For variety, we use the 
alternative discretization of (Vh) 2 .) We take the ends of the chain to be pinned; i.e. h (t) = 
= h N+ i(t). Again, the physics of the continuum problem admits only one asymptotic 
solution - that of a flat interface with hi(t) = 0. The explicit form of the discrete equations 
is 
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dhi/dt = -2/ii + h 2 + (X/A)h\ 

d/i i /dt = /i i _i-2/ ii + /i m + (A/4)(/i m -^_ 1 ) 2 , 2<i<iV-l (13) 
dh N /dt = h N -i - 2h N + (A/4)/i^_ 1 

The task of finding all the fixed points of the above equations is extremely difficult. We will 
content ourselves with studying an obvious fixed point: making the Ansatz hi{t) = h(t), we 
see that such a solution picks out two fixed points. The first is the physical one h(t) = 0, 
whilst the second is h(t) = 4/A. An analysis of the physical fixed point reveals, of course, 
its linear stability. The stability analysis of the unphysical fixed point follows from writing 
hi(t) = 4/A + hi(t). Keeping linear terms only, we have the stability equations 

dh l /dt = -2h x + 3h 2 (14) 
dh i /dt = h i - 1 -2h i + h i+1 , 2 < % < N - 1 (15) 
dh^/dt = 3/iAr_i — 2h,N (16) 

The stability criterion is now reduced to the following question. Does the matrix M defined 
by dh/dt = Mh have any positive eigenvalues, and if so, how do these eigenvalues depend 
on N. Details of a numerical solution are given in Appendix B. We find the following 
result: for iV < 6 there exists one positive eigenvalue, whereas for N > 6 there exist two 
positive eigenvalues which become degenerate and of order unity, in the limit of large N. 
We therefore see that this N body problem is also unstable and once again, the instability 
occurs for height deviations which become greater than a value of 0(1/ X). 

We shall leave a full discussion of these results for the final section. At this point we 
just reiterate the results of this section. We have studied the deterministic version of the 
conventional discretization scheme from a number of viewpoints. It has been found to 
yield both unphysical fixed points, and associated instability, whenever neighbouring height 
differences exceed a critical value of 0(1/ X). Thus, it has no connection to the continuum 
model in the strong coupling regime. The addition of noise to this scheme only lessens its 
range of applicability since (even for weak noise) strong fluctuations will create critical height 
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differences thus driving the evolution away from the correct asymptotic form. 



III. NEW DISCRETIZATION SCHEME 

The previous section has established the entirely negative result that conventional dis- 
cretization schemes of the KPZ equation are in fact physically unrelated to that equation, 
with the degree of physical integrity decreasing as one increases the coupling constant. The 
important question remains - how does one discretize the KPZ equation in such a way as 
to ensure stability and also to retain the important physics? In this section we provide an 
answer to this question which turns out to extremely simple, but also physically revealing. 

Before doing so, however, we consider the possibility that the discrete scheme displayed in 
Eq. (]4]) can be regularised by the addition of further terms, thereby avoiding any instability. 
Such an approach has recently been developed in the context of an exactly soluble mean-field 
model In its simplest form, it corresponds to the addition of a term /t(V/i) 2 V 2 /i to the 
right-hand side of Eq. ([[]), followed by a 'naive discretisation' as in Eq. (|4]). The resulting 
equation is soluble in the infinite-range limit where each height hi is coupled to every other 
height and, in the limit of weak noise, exhibits an interesting phase transition as a function 
of the variable g = A/ -J~kv . The strong-coupling phase g > g c is characterised by a double- 



peaked distribution for the fluctuations of the local heights hi from the mean height |T3 
This implies bumps of a definite height in the surface, reminiscent of numerical studies of 
a similar regularised model |14|]. This behaviour is, however, quite different from that of 
the supposedly related directed polymer model, leading one to question whether, even after 
regularisation, the discrete and continuum models belong to the same universality class. 
This is the question we now address. 

As we mentioned in the Introduction, there is a model intimately related to the KPZ 
equation; namely, the directed polymer in a random medium. The model representation 
consists of a linear diffusion equation with multiplicative noise, as shown in Eq.(|3]). The 
spatial discretization of this equation is completely under control, as the discretization only 
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affects the Laplacian term. The strong-coupling physics enters via the local noise term, 
and will not be affected by spatial discretization. The course is clear. We implement the 
spatial discretization via the diffusion- like equation and then make the inverse Hopf-Cole 
transformation on the lattice in order to obtain a discrete form of the KPZ equation which 
is guaranteed to be stable for arbitrary values of the coupling. 
So explicitly we rewrite Eq.@ as 

dwi/dt = iyja 2 ) ^{wj - w{) + (A/V) Wi rji. (17) 

j 

We now make the inverse transformation 

Wi(t) = exp[Xhi(t)/u] (18) 
which yields the new discrete form of the KPZ equation: 

dhi/dt = (z/ 2 /« 2 A)E{ ex P[( A /^)(^' - hi)) ~ + (19) 
j 

This equation appears unwieldy, but contains a very important message. It is useful to 
define local 'velocities' Vij = F(hi — hj) in such a way that one may write the discrete KPZ 
equation in the form 

dhi/dt = J2 v i,j+Vi- (20) 

j 

One then sees that the functional form of these velocities for the above discretization scheme 
takes the form 

F{x) = (is 2 /\a 2 ) [exp(-Xx/u) - 1]. (21) 

In contrast, the conventional scheme as shown in Eq.(|J) gives a functional relationship of 
the form 

F(x) = -(u/a 2 )x + (X/2a 2 )x 2 . (22) 

Comparing these two forms indicates that the conventional scheme is a truncated form 
of the new scheme to second order in (Xx/u) which will be accurate so long as x ~ (hi — 
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hj) < vj\ which is consistent with our findings in the previous section: the conventional 
scheme becomes unstable when neighbouring height differences occur which are of order 
vj\. The physics of the KPZ equation in the SC regime is such that large local maxima 
(i.e. those points which have large positive local height differences) are not evolved by 
the equation. They experience a weak negative force of order (v 2 /\a 2 ). We see that the 
conventional schemes fail because they actually boost such points in the positive direction. 
This unphysical behaviour will always be present if one attempts to discretize the KPZ 
equation using a finite number of terms in the height differences. The correct discretization 
schemes (of which Eq. fll9|) is one example) must be such as to freeze points which are local 
maxima. This effective freezing was already seen in the deterministic top hat problem 
studied in section 2. 

To end this section we make some practical remarks. The scheme presented in Eq. (|19D 
is unwieldly and will become more so when one discretizes time as well. In this case one 
must not discretize time in the h representation, but must inverse Hopf-Cole transform the 
space-time discretized equation for w. It is clear that one may as well numerically integrate 
the equation for w and then perform the statistical analysis on logarithms of that function. 
In fact one is essentially simulating the directed polymer problem. This implies that there is 
no easy way of directly integrating the KPZ equation. The only physically meaningful way of 
doing this is to use Eq . (|T9|) , which is tantamount to dealing with directed polymers. In fact, 
there is one example in the literature of a direct integration of the discretized equation for w 
in d = 2 ||15|| . The authors were motivated to integrate this form of the equation due to the 
reported problems in earlier work of numerical instability for larger values of the coupling 
constant. This motivation has been seen in the previous sections of the current work to be 
well justified. Their results for the dynamic exponents agree well with those found in ref. 
fT0|| , although they stress the difficulty of determining whether their results were obtained 
in the true asymptotic regime. It is certainly worthwhile to extend this type of numerical 
analysis to both larger values of the coupling constant, and to higher dimensions. 
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IV. NEW CONTINUUM ANALOGUES 



Aside from the practical concerns of how to integrate the KPZ equation, the previous 
sections have shed some light on what the SC physics of the KPZ equation actually means. 
What we have seen is that for large values of the coupling constant, local maxima in the 
interface become frozen or pinned. The interesting dynamical effects occur when we pump 
such a system with noise. Only the largest positive deviations will remain pinned for a long 
time, the smaller ones being swamped by the local noise fluctuations. In other words, the 
rare fluctuations become frozen in and therefore assume a much greater statistical weight. 

The question is how sturdy is this type of SC physics. Referring again to Eq.([H]), we 
see that the SC physics of the KPZ equation lies in the fact that F(x) is essentially zero for 
x > 0. The precise form of this function for x < may not be so important - this part of 
the function simply controls the rate at which lower points try to catch up with their frozen 
higher neighbours. To test the hypothesis that the KPZ universality class is determined 
only through the special features of F(x) for x > 0, it is interesting to consider other models 
which retain this structure. 

One such example is the continuum model defined via 

d t h = DV 2 h + r], (23) 

where 

I , V 2 /i< 

D = { (24) 

[Do, V 2 /i > 

This model is a not-too-distant relative of a stochastically pumped version of Barenblatt's 



equation |TT|,|T8 



The deterministic evolution of this equation consists of regions of negative curvature 
being strongly pinned, whilst regions of positive curvature are diffusively relaxed. In this 
sense it is an inherently SC realization of the KPZ equation. The relaxational time scales 
are certainly slower in this model than in the KPZ equation, but as we argued above, if one 
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believes that the SC evolution is dominated by frozen surface maxima, then this difference 
of time scales may not effect the universal features of the model. The pinning of regions 
of negative curvature in the above model actually makes it a prime continuum candidate 
for the description of the RSOS model, in which evolution is only possible within regions of 
positive curvature (although one must probably be more careful with the definition of the 
noise in this case.) It would be of interest to numerically study this model and ascertain 
whether it shares the same scaling features as the KPZ equation. Perhaps it is even possible 
to perform some analytic treatment. The nice feature of this model (as with many others 
which may be generated by altering the function F(x) for x < 0) is that it is intrinsically 
SC in the KPZ sense. We leave this speculative section now and end the paper with our 
more concrete conclusions. 



V. DISCUSSIONS AND CONCLUSIONS 



We have shown that conventional discretization schemes [p|4T0|| of the continuum KPZ 
equation (such as that shown in Eq.(f|)) have two failings which forbid them to remain faithful 
to the continuum equation, especially in the SC regime. Firstly they exhibit ghost fixed 
points (i.e. certain configurations which are stationary solutions) which are absent in the 
continuum; and secondly, these ghost fixed-points are unstable, so that the models deviate 
exponentially fast from the physical evolution. We believe we have clearly demonstrated 
these defects such that the previous results obtained from numerical integration of these 
schemes must be viewed as unsound. 

We have shown that by discretizing the directed polymer version of the KPZ equation, 
one may establish a stable, physically faithful scheme for integration of the KPZ equation. 
The main advantage of this new discretization is not so much its practical worth, but rather 
its illumination of the true nature of strong-coupling behaviour in the KPZ equation. On the 
basis of this new insight we have postulated a new criterion for models to be in the continuum 
KPZ universality class, and have suggested perhaps the simplest continuum member of this 
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class, which is related to the Barenblatt equation with a stochastic source. 

A perplexing question remains. Although we now understand how to make discrete mod- 
els faithful to the continuum KPZ equation, it is disconcerting that the discrete forms such 
as Eq.(f|) are not in the continuum KPZ universality class, when they seem so reasonable 
at first glance. It opens the possibility of a more tenuous link between microscopic (lattice) 
models and the continuum KPZ equation. What we have seen is that SC physics in the KPZ 
equation requires local maxima to be strongly pinned. This kind of physics is transparent in 
some microscopic models such as the RSOS model, and polynuclear growth models. How- 
ever, other microscopic models which might be better represented by naive lattice versions of 
the KPZ equation, will certainly not be in the true KPZ universality class. The cataloguing 
of the universality classes of the KPZ equation, and other non-equilibrium interface models, 
remains as ever an absorbing challenge for the future. 

TJN is grateful to Michael Swift for useful preliminary discussions, and acknowledges 
financial support from the Engineering and Physical Science Research Council. AJB thanks 
Matteo Marsili for discussions. 
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APPENDIX A: 



In this appendix we extend the results of the top hat problem studied in section 2. We 
shall connect the existence of a critical coupling constant to the known behaviour of the 
noisy KPZ problem in the vicinity of the weak coupling fixed point. First, we generalize 
the results to arbitrary dimension, by considering a pill-box in dimension d, of radius b and 
height H. Restricting our attention solely to the centre of the pill-box we have the exact 
solution 

h(0, t) = (u/X) In {1 + [exp(AV^) - 1] /(*)} (Al) 



where 



6/(4i4)V2 

m = TW) 1 d "" e ""'' <A2) 



o 

.2 



and T(z) is the Gamma function | 12fl . For vt ^> b this expression simplifies to 



v ' ; v 1 ' \ r(l + d/2) (47rz/t) d / 2 J v ; 

We see once again that there exists a critical value for the dimensionless parameter 
(XH/u). For (XH/u) <^ 1 the above expression may be expanded out to give the simple 
diffusion result for large times. Explicitly one has 

h{0,t) ~ b d H/{A-nut) d ' 2 . (A4) 

However, for (XH / 1/) > 1 we may only expand the logarithm for times larger than t* where 

t* = (b 2 /v) exp(2XH/du). (A5) 

For t 3> t*, the centre of the block once again diffuses away. However, for t <C i*, the centre 
is essentially frozen, and decays only logarithmically slowly. Explicitly one has 

h(0, t)~H- {dv/2X) \n(ut/b 2 ). (A6) 
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The message we wish to draw from this result is that for any finite dimension d, there 
exists a critical value of the coupling constant of order A c ~ (u/H), such that for A ^> A c , 
a surface deviation which has a scale of H will become essentially frozen, whilst a surface 
deviation of a smaller scale will diffuse away, almost unaware of the non-linearity If the 
reader feels uncomfortable with the discontinuous properties of the top hat function, we 
note that similar asymptotic results may be obtained by replacing the top hat over the 
range \x\ < b by an inverted parabola of central height H. 

Now consider the full noisy KPZ equation, starting from an initially flat interface. Over 
an initial time period the surface fluctuations will evolve according to the linear (or Edwards- 
Wilkinson IfLl^j) model. Accordingly the root-mean-square fluctuations grow in time as 



((Shy) 



2\l/2 



t^' A , d < 2 

(ln(t/t )) 1/2 , d = 2 (A7) 

const. , d > 2 

It is our contention that the SC physics of the KPZ equation will come into effect only when 
these rms fluctuations become of the same order as a critical height scale H, which is set by 
the condition \H/v ~ 0(1). This idea then leads to a prediction of the crossover time, t c 
from weak to strong coupling behaviour by equating ((<5/i) 2 ) 1//2 with H ~ u/X. One finds 
for d < 2 that t c ~ \~ 4 /( 2 - d ). For d — 2 one has t c ~ to exp[C(u 3 / DX 2 )] (where C is some 
number of order unity.) For d > 2 the rms fluctuations saturate to a constant value. One 
then has a very simple interpretation of the weak/strong coupling phase transition. For 
small A, the critical height H is large and thus exceeds the saturated rms fluctuations - 
hence there will be no possibility for the SC behaviour to set in. Conversely, for large A, the 
critical height H is small and will be quickly swamped by the weak rms fluctuations, such 
that the SC behaviour is initiated. 

The particular forms for the weak/strong coupling crossover time found above for d < 2 
have been found in previous studies by an explicit integration of the one-loop renor- 
malization group flow equations. The simple derivation above is useful as it brings out the 
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physical mechanism behind the flow to SC physics - i.e. the freezing of rms fluctuations of 
order vj\ and the subsequent non-linear evolution. 



APPENDIX B: 

In this appendix, we examine the eigenvalue spectrum of the matrix M that appears 
in the stability analysis of the A-site problem in section 2. Explicitly one has the A x A 
matrix 



M = 



V 



-2 3 

1 -2 1 

1 -2 














-2 10 

1 -2 1 
3-2 



(Bl) 



The interest lies in whether this matrix has positive eigenvalues (which signal an instability 
in the original problem), and if so, whether these eigenvalues survive in the large- A limit. 

We determine the eigenvalues by direct numerical computation using the Mathematica 
program. The result is that for increasing A, there exist two positive eigenvalues which 
become degenerate at a value of 0.1214(1) for large- A. A sample of the results for finite A 
is presented below: 



A 


4 


5 


6 


7 


8 


10 


13 


20 


25 


29 


h 


0.303 


0.236 


0.199 


0.175 


0.160 


0.141 


0.132 


0.122 


0.1215 


0.1214 


h 








0.000 


0.047 


0.091 


0.108 


0.121 


0.1212 


0.1213 
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Figure Captions 

Fig. 1: The time evolution of the centre of the top hat function, and of one edge of the 
top hat function, for A = 10A C . The curves represent the exact solution as given in Eq.(§), 
and the points represent numerical integration using the conventional discretization scheme 
shown in Eq.(Q). Parameter values are v = 1.0, a = 1.0, b = 10.0, H = 2.0 and for the 
numerical integration At = 0.002. For these values, A c = 0.5. 

Fig. 2: The phase space of the three-site problem. The dashed lines represent the domain of 
stability which has linear dimensions of 0(1/A). The insets indicate the site configurations 
corresponding to the various fixed points. 
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